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We further study the stochastic model discussed in Ref . |2| in which positive and negative 
particles diffuse in an asymmetric, CP invariant way on a ring. The positive particles hop 
clockwise, the negative counter-clockwise and oppositely-charged adjacent particles may 
swap positions. We extend the analysis of this model to the case when the densities of 
the charged particles are not the same. The mean-field equations describing the model 
are coupled nonlinear differential equations that we call the two-component Burgers 
equations. We find roundabout weak solutions of these equations. These solutions 
are used to describe the properties of the stationary states of the stochastic model. 
The values of the currents and of various two-point correlation functions obtained from 
Monte-Carlo simulations are compared with the mean-field results. Like in the case of 
equal densities, one finds a pure phase, a mixed phase and a disordered phase. 
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1. Introduction 



Phase transitions in stationary (far away from equilibrium) states present still many 
open questions. Some time ago we have introduced a simple one-dimensional three- 
state model which has interesting properties. In this model one takes a ring with L 
sites, the positive particles hop clockwise with a rate A: 

(+)(0)-(0)(+) (1.1) 
((0) is a vacancy). The negative particles hop counter-clockwise with the same rate A 

(o)(-)-(-)(o) (1-2) 

and oppositely charged adjacent particles may exchange positions 

(+)(-) -»• (-)(+) with rate g (1.3) 
(-)(+) -> (+)(-) with rate 1 (1.4) 

This model is translational invariant and the numbers of positive and negative particles 
are conserved. A detailed study of this model was done in Ref. |2| where we have 
considered the case in which the densities p (positive particles) and m (negative particles) 
were taken equal (p — m — p). Obviously 

p + m + v = 1 (1-5) 

where v is the density of vacancies. Using mean-field analysis, Monte-Carlo simulations 
and representations of the quadratic algebra given by the rates describing the processes 
( |l.l|H1.4j ) we have shown that the model has three phases. For q < 1 (independent on A 
and the densities) translational invariance is spontaneously broken (pure phase). In this 
phase, in the thermodynamical limit, the system organizes itself into a block of positive 
particles followed by a block of negative particles and one block of vacancies. All blocks 
are pinned in an arbitrary position on the ring. For 1 < q < q c (X, p) (called mixed 
phase) one has Bose-Einstein condensation in coordinate space (spatial condensation). 
In this phase one has again charge segregation but translational invariance is not broken. 
To understand better what happens (for finite but large L) it is interesting to start 
with the mean-field results. Neglecting fluctuations, the mean-field solutions, which 
break translational invariance, show two pinned adjacent blocks. One block (that we 
called condensate) contains only charged particles (no vacancies) with inhomogeneous 
distributions and a block (called fluid) which contains all the vacancies and where 
the particles have homogeneous distributions (see Ref. Fig. 25). As a result of 
fluctuations, the whole structure (condensate+fluid) makes a Brownian motion on the 
ring and the stationary probability distribution function is obtained taking with equal 
probability, the two adjacent blocks anywhere on the ring. For q > q c , one is in 
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the disordered phase; there is no charge segregation and the particles are distributed 
uniformly. In mean-field, one gets: 



function in the grand canonical ensemble as a function of fugacity a la Lee and Yang 
||, in order to get good estimates of q c (for other methods, see Ref. pf). 

The aim of this paper is to further study this model in the case of unequal densities. 
As we are going to see, the physics is different. Beside our physical motivation, we 
have also a mathematical one. The mean-field equations given by the present model, 
are coupled nonlinear differential equations and as we are going to show, we are able 
to find exact weak || solutions of these equations. We believe that without knowing 
the physics of the model it would have been very difficult to guess them. Since in 
the absence of vacancies, the coupled equations reduce to the Burgers equation 0, we 
call the mean-field equations "The two-component Burgers equations". In Ref. 0, for 
equal densities, we found stationary solutions of these equations. In the case of unequal 
densities we find the same three phases. Now however, the spatial structure (which 
exists in the pure and mixed phases) drifts with constant velocity around the circle. In 
the large L limit, the drift velocity vanishes exponentially in the pure phase, like 1/L 
for q = 1 and is finite for 1 < q < q c (\,p,m). 

From the mean-field results one can be tempted to conclude that one has no 
stationary solutions at all in the case of unequal densities even if we take into account 
fluctuations. Actually one can be left with this feeling from Ref. [0] where a different 
model was considered. That stationary solutions exist in our problem is known from the 
existence [||, § of a quadratic algebra with known representations and one can obtain 
in this way exact solutions for the probability distribution for the stationary state. 

In order to understand the physics of the problem it is instructive to have a look at 
the time evolution of a Monte-Carlo run which can be seen in Fig]]]. The parameters 
are chosen such that one is in the mixed phase. 

One notices that at a given time t one has two regions. In one of them (the 
condensate), there are no vacancies and that the boundaries are sharp. This allows 
a definition of the width of the condensate: one takes the last vacancy seen in the 
simulation on the left side and the first vacancy seen on the right side. In the condensate 
the particles are not uniformly distributed: the positive ones are predominantly on the 
left side and the negative ones are on the right side. The second region (the fluid) 
contains all the vacancies and the charged particles are distributed at random. If we 
know look at the figure downwards, one notices that the condensate drifts to the left 
with a constant velocity which can be estimated following for example the last vacancy 
on the left side of the condensate. Notice also that in the fluid, the positive particles 





As shown by one of us ||, one can use the distributions of the zeroes of the partition 
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Figure 1. A Monte-Carlo "film" for q = 1.2, A = l,p = 0.3, m = 0.1 (mixed phase), 
L = 100. The time t is running from the top to the bottom of the figure. The positive 
particles are depicted black, the negative ones-gray. 



drift to the right whereas the negative particles drift to the left. 

In Fig.||] we show the average values of the drift velocities obtained looking at the 
motion of the condensate as a function of the size of the system L in the pure phase 
(A = 1, q — 0.8), for q — 1 (A = 1) and in the mixed phase (A = l,q — 1.2). We have 
taken p = 0.4 and m = 0.1. 

In the pure phase, the velocity vanishes exponentially. A fit (also shown in Fig.0) 
to the data gives: 



V = 0.2742 exp(-0.02615L) 
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Figure 2. The drift velocity V as a function of the size of the system L. The upper 
curve corresponds to the mixed phase, the intermediate one to q — 1 and the lowest 
one to the pure phase. 

For q — 1, the drift velocity vanishes algebraically: 

V = 9.602/L (1.8) 
whereas in the mixed phase the fit gives: 

1/ = 0.116 + 6.9658/L (1.9) 

We now return to Fig.[l] and notice that the condensate does not only drift to the left 
but also does a Brownian motion. If one follows again the last vacancy on the left side 
of the condensate, it behaves like a random walker and therefore one can determine the 
diffusion coefficient D. This was done for different sizes of the system and the results 
are shown in Fig.^| for the same parameters as those used in Fig.|[ We notice that 
similar to the drift velocity V, for large values of L, the diffusion coefficient converges to 
a constant value in the mixed phase, vanishes algebraically for q = 1 and exponentially 
in the pure phase. 

As a result of these observations, it is clear that in the pure phase, in the 
thermodynamical limit, one gets three adjacent pinned blocks and like in the case of 
equal densities, translational invariance is broken. The situation is different in the 
mixed phase. If we follow again the last vacancy on the left side of the condensate, 
since it moves like a random walker with a finite diffusion coefficient, the dispersion 
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Figure 3. The diffusion coefficient D as a function of the size of the system L. The 
upper curve corresponds to the mixed phase, the intermediate one for q — 1 and the 
lowest one to the pure phase. 

increases with time up to when it reaches the perimeter of the ring. For later times one 
can find this vacancy anywhere on the ring. This implies that in the stationary state, 
one finds with equal probability the structure condensate+fluid in any position on the 
ring. The existence of such a structure obviously implies charge segregation (long range 
correlations) without breaking of translational invariance. 

We can understand better now what is going to be the role of mean-field calculations. 
On one hand, one can estimate non-equilibrium quantities like the drift velocity, on the 
other hand one can obtain the condensate+fluid structure. Once the latter is known, 
one can obtain the probability distribution of the stationary state taking as mentioned, 
with equal probability the charge distribution obtained in mean-field in any position on 
the ring. 

The paper is organized as follows. In Sec.|2] we consider the mean-field equations 
for the case q > 1. Those are the two-component Burgers equations. In this section 
we shortly remind the reader some of the results obtained in Ref. for the case of 
equal densities and also present some new ones. In Sec.| we find exact weak solutions 
of the two-component Burgers equations for the case of unequal densities which tell 
us under which conditions the condensation phenomenon appears. In this way one can 
obtain the function q c (X,p, m). Moreover one obtains the drift velocity V and the charge 

5 



- 






















- 


- 

- 






















- 

- 






x 


























X 

Tr x 






• X X X. 


....X X.... 


....X.... 


....X X 








k 




















X. i 


~ + 






























X 


-■x-- 






















X. , 1 , 1 


i 


— X— 


-x— - X - - X - 

1 , 1 






— x- 




■-Xr i 



distributions in the condensate. In Sec.f| we discuss the case q = 1, this case is much 
easier to study than the q > 1 case. We have not considered the case q < 1 since it can 
be studied along similar lines (see Ref. for the case p = m). 

In Sec.[5] we discuss the results of Monte-Carlo simulations and compare them with 
the results obtained in mean-field used as described above. We also give some results 
using the algebraic approach |2| which is, in the case of unequal densities, harder to 
apply. The conclusion is presented in Sec.|6|. 

Since this work is the continuation of the one presented in Ref. 0, this paper is 
not self-contained and it can't be understood without first reading Ref. unless one 
is interested only in the solutions of the two-component Burgers equations described in 
Secs.|, | and |. 

We would like to mention that the model described above was subsequently studied 
again on a ring, in the case in which only one vacancy (seen as a defect) is present [[/J 
and for a special case of open boundaries || . 



2. The two-component Burgers equations (q > 1) 

The mean-field equations in the continuum of the model defined by Eqs.( |l.l|) -( |1.4|) are 
(Ref. [| Appendix B): 

where p± are the local densities of the charged particles (po is the local density of the 
vacancies) satisfying the obvious relation: 

p+ + p_ + po = 1 (2.2) 

j + is the local current density of positive particles, j_ is minus the local current density 
of negative particles (the minus sign is convenient if one looks at stationary states in the 
case of equal densities since with this choice of signs one gets j + = jJ). The expressions 
of the current densities are: 



j± = Ap ± p + (g-l)p+p- + ^ 



/ _\_ i\ f d P- d P+\ i\( d Po d P± 

(<? + !) P+^r ~ p-^t ) ±x \ p±^z - Po- 



(2.3) 



dx dx ) \ dx dx 

We are looking at periodic solutions of the Eqs. (|2.1|): 

P±(x,t) =p±(x + L,t) (2.4) 

where L is the perimeter of the ring. We will also be interested in solutions with a drift 
velocity V counter clockwise (see Fig.|l|): 



p±(x,t) = p±(x + Vt) 
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(2.5) 



which implies that in the moving frame, the quantities 

C± = ±Vp ± + J± (2.6) 
are independent of x. It is useful to define 

V=2 0<y<l (2.7) 

and 

J± = 7 , W = , l± = , v = — , 77 = , p = — 2.8 

q — 1 q — 1 g — 1 2L(g — 1) g — 1 2L 

In Eq.(|2.8|) we have assumed that q > 1. Rescaling the time, instead of Eqs. (|2.1|),(|273|)- 
( f2.6| ) one finds: 

«± = W±Po + ^P+P- + f P+^ P-^~ =•= P P±^ ^"a - 2 - 10 

V dy dy J \ dy dy J 

t±(y,t)=£±(y+l,t) (2.11) 

e±(y,t) = t±(y + wt) (2.12) 

J± = ±Wp±±£± (2.13) 

where u; = 1 . Since the numbers of positive and negative particles are conserved, we are 
looking for solutions for given densities p and m: 

1 rl 

p+ dy = p, / p_dy = m (2.14) 
jo 

We call the Eqs. ( J2.10|) — the two-component Burgers equations. We have coined 
this name to the equations since if the density of vacancies po is equal to zero, we can 
make in Eqs.( |2.lOD the substitution 

*±=l±f ( 2 - 15 ) 

and see that / satisfies the Burgers equation || 

df d 2 f df 
dt dy 2 dy 

where ( = 2 and v is the damping constant whose physical interpretation in our model 
can be read of from Eq.( |2.8j ). One notices that in the infinite volume limit, the damping 
constant vanishes and one obtains the inviscid Burgers equation ||. The constant p 
appearing in Eqs. fl2.lOD will be called moisture constant. Notice (see Eq.(|2.8|)) that the 
moisture constant also vanishes in the large L limit. 
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An important result of our investigation is that one is able to give exact solutions 
to the two- component Burgers equation. Let us review the results obtained in Ref. [[J 
(Section 6). For equal densities, in the mixed phase, inspired by the results obtained 
numerically (see Fig.25), one can obtain the stationary solutions in the following way. 
We divide the segment < y < 1 into two smaller segments: < y < a where we have a 
condensate (no vacancies) with a non-uniform distributions of the charged particles and 
a second segment of length b, (a + b = 1) where the particles have a uniform distribution 
(fluid). Therefore in the condensate one can use the Burgers equation whereas in the 
fluid one has almost nothing to compute. Numerical investigation teach us how to sew 
the two segments: the density of positive particles is smooth at the right-side of the 
condensate (y = a) but not at the left-side (y = 0) where it is discontinuous. For the 
negative particles the discontinuity is at y = a as a consequence of the CP invariance of 
the problem. The solutions obtained in this way are weak solutions of the differential 
equations which correspond to taking the moisture constant p equal to zero. They 
represent an excellent approximation not only of the concentration profiles obtained 
numerically for the mean-field equations but also of the concentration profiles seen in 
computer simulations (see Fig. 8 and 9 in Ref. 0). 

The values of a are dependent on u, rj and the density p = m = p. The domain in 
the space of the variables u, rj and p where solutions exist define the mixed phase. In 
the thermodynamic limit [y = 0) one obtains for b the value 

B = t + 2 "' (2-17) 
2(1 + 'f) ^ ' 



The result mentioned in section |I]for q c (see Eq. ( |1.6| )) is obtained from Eq. fl2.17|) when 



takes B = 1 (the condensates disappears). To the results obtained in Ref. [Q] we would 
like to add some supplementary data which are relevant also for the case of unequal 
densities. For a given value of rj there is a maximum value of v (z/ max ) beyond which 
the condensation does not take place anymore. Moreover, one would expect that for 
^max the length of the condensate a is zero and that a increases as the damping constant 
decreases. This is not the case: for v = z/ max one gets a non-vanishing value for a equal 

to O m i n . 

In FigsfJ and ^ one can see this phenomenon. We work with the variables q, A 
and L instead of v and r] (three variables instead of two) not only because our physical 
intuition is in these variables but also because in section |5| we will look into the stochastic 
problem where, unlike in mean-field with /i = 0, on has all the three variables. 

One notices that as we approach the critical point (q c = 1.5714 as obtained from 
Eq.(|1.6|)), the minimal length L min diverges and the the smallest condensate a min 
vanishes. 

We should stress again that the solutions described above are weak solutions for the 
stationary states of the two-component Burgers equation. The current is constant on 
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Figure 4. The minimal length of the system to have a condensate for different q and 
A = l,p = m = 0.2. The data is very good fitted by L = I7(q c — q)~ 3 / 2 . 

the segment < y < 1 but the densities are discontinuous. One can ask the question 
how do these discontinuous solutions look like if one looks at the "microscopic" problem 
on the lattice both in mean field and in Monte-Carlo simulations. As one can see in 
Ref. H, Fig.25 or in the present paper Figs. and 0, there are two dramatic changes 
in the densities over not more than a few sites for any lattice size L. If one works with 
the y = x/L variable, these localized, abrupt, changes lead to discontinuous solutions. 

One can ask if in the mixed phase the condensate+fluid solutions are the single 
ones. The answer is no. Obviously one has a solution where all the densities are 
constant (pure fluid). More interestingly, for a given value of u, 77 and p one can have 
solutions with n condensates (fluid+condensate-h ■ -+fluid+condensate) all of them with 
the same distribution of charged particles and located anywhere on the ring. It is easy 
to show that in order to get the solutions with n condensates one has to use the solutions 
obtained for a single condensate in which we take the damping constant equal to nv. We 
didn't study the problem of the stability of these solutions. What we do know is that 
we have also obtained numerically these solutions by using the time-dependent mean 
field equations (see |2| Eq.(B2)) and that we have obtained, depending on the initial 
conditions, solutions with no condensate, with one condensate and with two condensates. 
We didn't take large lattices where one could have seen more than two condensates. 
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Figure 5. The length of the condensate when condensation starts for different q and 
A = l,p =m = 0.2. 

One comment is in place concerning the multi-condensate solutions of mean-field. In 
the presence of fluctuations, each condensate behaves like a random walker who sooner 
or later will meet another random walker. When two condensates meet, they melt into 
a single condensate so that if one starts with n condensates, we end up with a single 
one. This is seen in Monte-Carlo simulations. 

3. Roundabout solutions of the two-component Burgers equations (q > 1) 

Before describing roundabout solutions of the differential equations (|2.9| ) obtained in 
the case of unequal densities for the positive and negative particles, it is interesting to 
give some results obtained looking at large time solutions of the mean-field equations 
on the lattice (Ref. [0 Eq.(B2)) obtained numerically. In Fig. 6 we plot the current 
densities as a function of the densities of the positive and negative particles for a lattice 
with 1000 sites for q = 1.2, A = 1 and p = 0.4 and m = 0.1. We are going to see later 
that we are in the mixed phase. We notice that in agreement with Eq.( |2.7| ) one obtains 
straight lines and from their slopes one can with high precision the drift velocity V. 

In Fig.0 we show the distributions of the current densities and of the particle 
densities for a shorter lattice (L = 200) keeping all the other data unchanged. Like 
in the case of equal densities, one observes a condensate (a region without vacancies) 
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Figure 6. The current densities as a function of the densities for q = 1.2, A = 1. 
L = 1000, p = 0.4, m = 0.1 (lattice mean-field results). Using Eq.( ^.6[ ) one finds 
V = 0.11878. 



and a fluid with uniformly distributed charged particles. Based on these numerical 
observations one can start to look for roundabout solutions of Eqs. ( |2.9| ), taking into 

account Eqs. (|2.11|) -( p.l4|) . We first start with the condensate (0 < y < a) where we 
have no vacancies. Denoting 

1 + W 

pl = —^±u (3.1) 

J = J -^- (3.2) 

from Eq. fl2.13f ) we obtain: 

J + -J_ = W (3.3) 
and from Eqs.( |2.13"D and Q2.10| ) we get: 

a 2 =A J_ 1 _ W 2 = _ 4(m 2 + ^ (34) 

dy 

We have denoted by p± the densities of the charged particles in the condensate. 
Obviously 

P+ + P- = l (3-5) 
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Figure 7. The profiles of the densities and currents (scaled by a factor 3) obtained 
numerically on a lattice with 200 sites. The other data are like in Fig|| The condensate 
is drifting to the left. 



We can integrate Eq. (|3.4|) and find: 

u = /3u ta.n[P (aip — y)] (3.6) 
where (p is an arbitrary constant and 

P = f (3-7) 

At this point J, W and if are unknown constants which still have to be fixed. We 
now consider the fluid (a < y < 1) where all the vacancies are concentrated. In this 
domain the densities p± and p{j are constant (independent of y) and verify the relation: 

P% + P*+pl = l (3.8) 

We obviously have 

Po = v/b (3.9) 



where v is the density of vacancies and b the length of the fluid. From Eq.( 2.10 ) we get: 
£l = VPlPo+PlP- (3-10) 

We denote 

A = p R + -p R _ (3.11) 
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and get 



and 



W = V Ap R 

a 2 = -A 2 (r? - l) 2 - (p*) 2 (2r? - 1) + 2p*( V - 1) 



(3.12) 



(3.13) 



We now sew the condensate with the fluid asking for the following conditions to be 
fulfilled: 



P c + (a)=pl, p1(0)=P- 
Using Eqs.O and we get: 



, n , 1 w fl , . 1 w fl 
w(o) = g - y - p- » = 2 + Y - p+ 



and therefore: 



pfl = [3v( tan[/?a(l - <p)] + tan[/3a<p]) 



and 



- A = /3z/(tan[/3a(l - 9?)] - tan[/3a<p]) 
We now use Eq.( |2.14D and get for n condensates: 

pa pa 

n / p c + dy + p+b = p , n / p^Ldy + p®_b = m 
Jo Jo 

It is useful to denote 

d = p — in 
and from Eq. ( |3.18| ) we obtain: 

d=(Wa + 2f udy)n + Ab 
Jo 

We have to make again a change of notations: 

7T 7T 

P a( p = 2 - ^1 , M 1 - V) = 2 ~ ^ 2 

In the new notations Eqs. (|3.16|) , (|3.17|) , ( |3.20|) and ( |3.13|) become: 
Po = i 71 ~ + ^2)) -[cot ^2 + cotV'i] 

(X 

A(r]-1) = (tt- (^i + ^ 2 ))-[cot^ 2 -cot^i] 



d — A(r/ — (77 — 1)6) = 2i/n log 
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/ sin ip 2 
\ sin ipi 



(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 
(3.23) 

(3.24) 




Au 2 

— (vr - fa + = -A 2 ( V - l) 2 - {plf(2 V - 1) + 2rffo - 1) (3.25) 

CL 



Taking into account Eq.( |3.9| ), the four unknown b, ipi, ip2 and A can be obtained from 
the four equations fl3.2243~25l) once p and m are given. In this way one can determine the 



profiles of densities in the condensate and the fluid, the currents and the drift velocities 
once we give 77, u, p and m . We have not determined the domain of these variables 
where the solutions exist. Having in mind the stochastic process where we have used the 
mean-field solutions in order to explain the properties of the stationary distributions, we 
have just looked at the values q = 1.2, A = 1, p = 0, 4 and m = 0.1 in the mixed phase. 
As mentioned already before, and as can be seen from the equations given above, if we 
know the solutions for one condensate for any value L, one can obtain the solutions of 
the problem for n condensates for a lattice size L taking the solution for one condensate 
for a lattice of size L/n. 

In Fig.|8| we show the drift velocities for the condensates which appear up to 
L = 1000. One notices first that the first condensate appears for L = 60. This implies 
that that one obtains n condensates when (roughly) L > 60n. Next one can see that 
for a given value of L many condensates move faster than fewer ones. 

Various applications to the stochastic process of the distributions of the charged 
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particles obtained solving the transcendental equations given above can be found in 
Sec.|5]. (see Figures and Pp . 

The calculations simplify in the large L limit when the friction constant v constant 
is small. Denoting the limiting values by 



and 



we obtain 



B, A^D 

C = (77-1)5 
270/(77 - 1) 







and 



D 



1 - B 



d 



Pa 



— — D 

C 



v 
B 



^2 



2m/(r] - 1) 
1-5 



rj-C 

C can be determined from the cubic equation 



C 3 -C 2 



v , d 2 

2n + -(77 - 1) H 

' 2 W ; 2v 



2 

77 7J 



+ C77 [77 + 7; (277-I)] - -^—(277 



r 







(3.26) 
(3.27) 
(3.28) 
(3.29) 

(3.30) 



We have used the equations (|3.29|) and ( |3.30|) in the case q = 1.2, A = 1, p — 0.4 and 
m = 0.1 and have obtained 



V = 0.11180 , B = 0.57929 , D = 0.11181 
0.86313 



Pi 



j% = 0.1076 , j fl = 0.0111 , jl = f_ = 0.03437 



(3.31) 



Notice that the value obtained this way for V is closed to the value determined 
numerically on the lattice as given in the caption of Fig.|| As we are going to explain 
in Section |5], the average values of the current densities over the ring are very useful. In 
the large L limit they are given by 

( J± ) = (l-B)f ± + Bjl (3.32) 
Using the values given in ( |3.31|) , one obtains: 

( J+ ) = 0.07678 , = 0.0288 (3.33) 
These values are going to be used in Section |5|. 



4. Mean-field solutions for q = 1 

The starting point is again the Eqs.( [2.5|) and ( |2.6| ). Making the notations 



V= L> 



P 



2XL 



7+ V 
— W = — 
A A ' 



/ - C± 



(4.1) 
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we find the expression (|2.10| ) with r\ = 1 and to = 0. We repeat the procedure used in 
the last section assuming again that the moisture constant fi is zero. This assumption is 
again based on numerical solutions of the mean-field equations as can be seen for example 
in Fig. 9 of Ref. in the case of equal densities. In the region without vacancies, using 
again the substitution ( |2.15| ) we obtain the Burgers equation ( |2.16| ) with ( = 0. 
For the condensate we use again Eqs.(|3.1|)-( |3.3| ) and get: 



J 



J+ + J- 



v- 



dp c . 



(4.2) 



From Eq. 

(0 < y < a) 

where 
and 



2 2 ' ' r+ " dy 

we obtain the distribution of the positive particles in the condensate 



A Wy 

Ae~ + J 4 



.4 



-p+(0) +p+(a) 

Wa 

e " — 1 



W 



e - a p+(0) - p + (a) 



Wa 

e " 



We also have: 



J = W 



y + (0)-p c + (q) 

Wa 

e v — 1 



In the fluid (a < y < 1) where the densities are constant, we obtain 
J + - J = W(l - pS) + (pi - p fl )pg 



and 



as well as: 



J 



pDpI 



pi - P R =w 

Since J has to be constant over all the segment < y < 1, we get from Eqs.i 
& 



W 2 + (l- P l) P l = -W + 2Wp\ - 2u 



dp% 
dy 



(4.3) 
(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 
and 

(4.10) 



We now ask that the density of positive particles is smooth at the right hand side of 
the condensate (y = a) and that the density of negative particles is smooth at the left 
hand side of the condensate (y = 0) and get 



p c Ja) 



-(1 + W 



Pi) 



4(0) = -(l + W + p5) 



(4.11) 
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as well as: 



J 



<-> ' r> Wa PO 

z z e v — l 



Wa 

fl e " + 1 



e " 



and 



1 



4(v) = 5(1 + wo + + *W 

z 2 e f — 1) 



(4.12) 
(4.13) 

(4.14) 



We now use Eq. ( |3.18|) which fixes the densities of the charged particles and use Eq. (|4.14j ) 
to obtain: 



PoC 1 -a)=v 
p - m = W + Pq 



Wa 

e~ + 1 

t Wa ~ 

e « — 1 



2i/ 



(t> is the density of vacancies). It is convenient to denote 
Wa 



v 



X 



From Eqs. (|4.13|) and ( |4.16 ) we obtain two equations for the unknown x and a: 



p — m 
1 



xv va 

— + 1 

a 1 — a 



e x + 1 2 
e x — 1 x 



1 — a 



xi/ + 1 
a e x — 1 



(4.15) 
(4.16) 

(4.17) 

(4.18) 
(4.19) 



[we have taken into account the relation ( |4.15| )). 
In the large L limit (very small v) we obtain 



Pi 



(4.20) 



which means that there are no charged particles in the fluid therefore the length of the 
fluid is 



B = v 

From Eq. Q4.18] ) we get 

p — m e x + 1 2 
p + m e x — 1 x 

which gives x once the densities p and m are given. Knowing x we get W: 

xv 



W 



p + m 



(4.21) 



(4.22) 



(4.23) 
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Figure 9. and 3(j_) from Monte-carlo data for q = 1.2, A = l,p = 0.4, m = 0.1. 

the lines come from mean-field in the L — > oo limit. 

It is an interesting exercise to derive the drift velocity V from these equations. For 
p = 0.4 and m = 0.1 we get V = 9.602/L in excellent agreement with the the estimate 
obtained from Monte-Carlo simulations (|1.8|) . One can easily derive the average (over 
the ring) values of the currents: 

(j+) = ^—(^-^-p) , {]-) = —- ^-T +m 4 ' 24 
p + m \e x — 1 / p + m Ve 1 - 1 / 

and notice that they decrease like the inverse of the length L of the chain. 

5. The mixed phase. Comparison between mean-field results and 
Monte-Carlo simulations 

In the last two sections we have shown, the mean-field equations related to the stochastic 
process defined by Eqs.( |l.l| )-( |1.4j ). In Ref. we have discussed in detail, in the case of 
equal densities, how to use mean-field in order to explain the properties of the stationary 
states. One assumes that in the mixed phase, one can take with equal probability in 
any position of the ring the structure (condensate+fluid) obtained in mean-field. In this 
section we extend this discussion for the case of unequal densities. We limit ourselves 
to the case p = 0.4, m = 0.1, and A = 1 only. 

In Fig.|S] we give the average values of the currents and (J-), for q = 1.2 and 
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Figure 10. The current of positive (top), negative particles (bottom) from Monte- 
Carlo Simulations and from mean-field calculations in the limit L — > oo (dashed lines) 
and homogeneous mean-field (solid lines). The parameters are A = 1, p = 0.4, m = 0.1 
and L = 800. 



various lattice sizes. Within the errors they look like they converge to the large L limit 
obtained in mean-filed. Since in mean-field the current-densities are not constant (see 
Fig.0), one has to compute the average value of the current densities j + and j_ on the 
ring. For the case q = 1 the expressions of these averages are given in Eq. (|4.24|) . For 
q > 1, in the infinite L limit only, the calculation of the averages is again simple since 
in the condensate the positive and negative particles have a uniform distribution like in 
the fluid (see also Eqs. (|3.32|) -( |3.33|) ). The calculation is more complicated in the case 
of large but finite values of L. 

The average values of the currents for various values of q are shown in Fig. 10. The 
mean-field values are obtained as described above, the Monte-Carlo values are obtained 
using results from finite lattices and extrapolating to L infinity. Similar to the case of 
equal densities, the currents are very well given by mean-field up to a value q c (around 
1.6). This is the mixed phase. For q > q c the currents have a different behaviour (one 
is in the disordered phase). 

In Figs.|TT] and [T^ we show the two-point correlation functions as a functions of the 
distance R = yL obtained in Monte-Carlo simulations (L = 200) and in mean-field. In 
the case of correlation functions, in order to obtain the mean-field results one had to 
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Figure 11. The correlation functions co,o (top) and c+,_ (bottom) for q = 1.2, A = 
1, p = 0.4, m = 0.1 and L = 200. Each curve consists of the Monte-Carlo and mean- 
field data lying perfectly upon each other. 



use the results of Sec.|] for L = 200. The agreement between the simulations data and 
the mean-field results is remarkable. 

Since one is in a stationary state, one can't talk about drift velocities which are 
seen in non-equilibrium or in the roundabout solutions of the mean-field equations. 
Let us observe however that when we use the mean-field results in order to derive the 
correlation functions given in Fig.0, for small values of y, one obtains the following 
value of the drift velocity V 



V = A lim 



c +,o — c o,- 



j/->0 1 — p 



m 



(5.1) 



This expression is so simple because in this limit, only the fluid and not the condensate 
plays a role. Using the results of the Monte-Carlo simulations and Eq. ( |5.1| ) one can 
estimate the drift velocity. One gets V = 0.14 which is perfectly consistent with the 
value obtained looking at the random movement of the condensate (see Fig.|^ and Fig.|l~3|) . 

Up to now we have discussed only the Monte-Carlo data. In Ref. we have 
also shown how to use the algebraic approach in order to obtain the current using the 
grand canonical ensemble (for correlation functions this approach gives wrong results). 
In the case of unequal densities one has to deal with two chemical potentials and the 
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Figure 12. The correlation functions c+ o (top) and Co,- (bottom) for q = 
1.2, A = 1, p = 0.4, m = 0.1 and L = 200. Each curve consists of the Monte- 
Carlo and mean-field data lying perfectly upon each other. The drift velocity is 
A limj,^o( c +.o — c o,-)/(l — P — fri)- For equal densities all four curves are lying on 
top of each other. 



calculations become difficult. In Fig.[L3]we give the values for the quantity 

V = {3+) - {3 - } (5.2) 
1 — p — m 

obtained using the algebraic approach and compare it to the corresponding values 
obtained in mean- field. In mean- field V as defined by Eq.( p.2| ) is indeed the drift velocity 
(use Eq.([2.6|) and compute C± in the fluid). The agreement between the two calculations 
looks reasonable. 



6. Conclusion 

In this sequel of Ref. we have introduced what we call the two-component Burgers 



equations defined by Eqs.( |2.9| ) and (|2.10| ). These equations correspond to the mean- 
field approximation of the stochastic model defined in Sec.|l] We have found roundabout 
weak solutions of these equations for d = 1 in Sec.^ and for rj = 1 and u = in 
Sec.f|. These two cases correspond to the mixed phase, respectively to the q = 1 case, of 
the stochastic model. The weak solutions of the two-component Burgers equations are 
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Figure 13. The "drift velocity" for q = 1.2, A = 1, p = 0.4, m = 0.1 from the 
alegbraic approach V = — (j_))/(l —p — m) (points) and mean-field (line). The 
straight line gives the mean- field value for the limit L — > oo. 



obtained if we take the moisture constant /x = in these equations. When compared 
with the solutions of the equations when the moisture constant is not zero, which can 
be obtained numerically, one sees that one gets a very good approximation. 

The analytical roundabout solutions presented in Sec.||] allow to determine the 
domain in the parameter space where the mixed phase exists. Moreover, using these 
solutions one can obtain predictions for the stochastic problem as explained in Sec||. 
These predictions are in good agreement with the results of Monte-Carlo simulations on 
large lattices, also for the case in which the densities of positive and negative particles 
are not equal. Since Monte-Carlo simulations are much time-consuming, in the present 
paper we have performed a less careful analysis of the data as compared with the case 
of equal densities. 

Before finishing this paper we will like to comment on a paper |H| where, in the 
case of equal densities, one questions the existence of the mixed phase. If one uses the 
grand canonical ensemble and the algebra presented in Ref. || and two hypotheses one 
argues that for " cosmological" large lattices the mixed phase disappear. For example, 
if one takes A = 1 and q = 1.11, the dimension of the lattice should be of the order 
10 70 in order to "loose" the mixed phase. If the value of q gets closer to 1 one get 
lattice lengths of the order 10 490 or more. If the calculations of Ref. ^TJ describe indeed 
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the reality this would be a fascinating phenomenon, certainly more interesting that the 



model itself. In Ref. [10| one has made one mathematical hypothesis which was proven 
correct [IT] and a technical hypothesis according to which one neglects in the calculation 
configurations with no vacancies. This second hypothesis, which was used already in 
Ref. 0, is probably also irrelevant |jT2"|| . Our own doubt about the results of Ref. [PHI 
come from the use of the grand canonical ensemble. As stressed in Ref. 0, for stochastic 
processes as defined for example in Sec.|T|, there is no grand canonical ensemble. One 
can make an ad hoc definition, hope to be lucky and get the right results. As shown in 
Ref. ||, for the two-point correlation functions in the mixed phase, the canonical and 
grand canonical ensembles give very different results (see Fig.21). It looked however 
that for the one-point function (the current) the results were correct. We don't want 
here to defend here the existence of the mixed phase (if it survives up to lattice lengths 
of the order of 10 70 , it is enough for physical purposes), on the contrary, we think that 
the real challenge is to clarify the problem. One more argument in the favour of the 
existence of the mixed phase can be found in Ref. |13| where one considers the open 
system (no grand canonical calculation in this case) for the same bulk parameters as 
those taken on the ring when one has the mixed phase. The input and output rates are 
chosen symmetric for the positive and negative particles. One sees a first-order phase 
transition between a phase where the density corresponds to the fluid (in the mixed 
phase on the ring) and a maximum current phase which corresponds to the condensate 
on the ring. (We remind the reader that the densities in the fluid depend on q and A 
Only and not on the density p = m = p). 
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